CMB and matter power spectra from cross correlations 
of primordial curvature and magnetic fields 



Kerstin E. Kunze0 

Departamento de Fisica Fundamental and IUFFyM, 
Universidad de Salamanca, 
Plaza de la Merced s/n, 37008 Salamanca, Spain 



Abstract 

A complete numerical calculation of the temperature anisotropics and the polarization of the cosmic 
microwave background (CMB) is presented for a non zero cross correlation of a stochastic magnetic 
field with the primordial curvature perturbation. Such a cross correlation results, for example, if 
the magnetic field is generated during inflation by coupling electrodynamics to a scalar field which 
is identified with the curvaton. For a nearly scale invariant magnetic field of 1 nG it is found that 
at low multipoles the contribution due to the cross correlation dominates over that of the pure 
magnetic mode. A similar behaviour on large scales is found for the linear matter power spectrum. 



1 Introduction 

The presence of magnetic fields on very different scales has been confirmed by observations [lj . On 
galactic scales there is an abundance of magnetic field detections with typical field strengths in the 
fj,G range at present time. Over recent years there have been claims of truly cosmological magnetic 
fields, that are not associated with any virialized structures such as galaxies or clusters of galaxies 
[2]. However, this has been challenged in [3]. 

There are many models to explain the origin of large scale magnetic fields, for reviews see, 
e.g., [4]. In particular, in the early universe conditions could have been such as to naturally 
generate magnetic fields strong enough to seed a galactic dynamo. The subsequent amplification 
since galaxy formation then results in [iG magnetic fields today. Inflationary models seem to be 
especially attractive to generate the initial seed magnetic field since the correlation lengths can be 
very large. The problem usually is with the magnetic field strength which in spatially flat models 
is too small to account for the initial magnetic seed field due to the global conformal invariance of 
standard electrodynamics in these backgrounds. However, in spatially curved models this is not 
the case and thus not a problem [5] (see however, [6]). 

Therefore, when considering the standard flat ACDM model conformal invariance of electro- 
dynamics has to be broken which in the simplest case is realized by coupling to a scalar field [7]. 
However, there are possible problems with back reaction effects onto inflation pHE)]. Other possibil- 
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ities include coupling electrodynamics to curvature [TU] of which some models also have problems 
due to ghost instabilities [11] . There are models in which magnetic fields of cosmologically relevant 
strengths can be generated during inflation |10} I12j . However, since in this case inflation is of 
power-law type there is yet another constraint coming from the fact the curvature perturbation 
for the relevant parameters cannot be created during inflation. It is possible to use the simplest 
curvaton scenario to generate the curvature perturbation after inflation [13]. This, however, results 
in a further restriction of the parameter space but still allows to generate magnetic fields strong 
enough to act as seed fields for the galactic dynamo [13] . 

In [15] a particular model has been studied in which magnetic fields are generated during 
inflation by coupling electrodynamics to a scalar field which is not the inflaton. In this case it has 
been shown that there is a non trivial correlation between the scalar field perturbation and the 
magnetic field [15]. Moreover, it was pointed out that identifying the (spectator) scalar field with 
the curvaton determines the cross correlation of the magnetic field with the curvature perturbation. 
This is the case to be considered here. The aim is to calculate the anisotropies of the temperature 
and polarization of the cosmic microwave background (CMB) due to the cross correlation between 
the magnetic field and the primordial curvature perturbation. Upto now, in general, the effect 
of magnetic fields present before decoupling on the CMB have been calculated assuming either 
complete correlation |16j-|19j or no correlation at all \20\ [T9] between the primordial curvature 
perturbation and the magnetic field. In most works the magnetic field is assumed to be non helical. 
However, the helical case has also been studied which is particularly interesting since it induces odd 
parity correlations of the CMB modes |21] . In |22] the effect of a non vanishing cross correlation 
between the curvature and the magnetic field due the evolution of the magnetic field has been 
studied. In that case the effect on the CMB temperature anisotropies and polarization has been 
found to be much below that of the compensated mode for a magnetic field which only redshifts 
with the expansion of the universe. 

The magnetic field does not enter linearly in the perturbation equations but rather in form of 
its energy density and anisotropic stress which are quadratic in the magnetic field implying that 
non gaussianity is induced even at linear order. The resulting bispectra have been calculated [23] , 
Thereby the cross correlations between the curvature mode and the magnetic field contributions 
are determined by the bispectra of the magnetic field and the curvature mode. In [24J the bispectra 
induced by coupling of the inflaton to electromagnetism have been calculated. This model has been 
also considered in [25] where the corresponding CMB anisotropies have been calculated. 

2 The model 

The background metric is assumed to be of the form 

ds 2 = a 2 (ri) (-dif + Sijdx^x^ , (2.1) 

where a{rj) is the corresponding scale factor. 

The model of magnetic field generation during de Sitter inflation used in [15] which is based on 
[El [26] is described by the action 

S = J d 4 xV^ (—W^F^F^ - \ {d4>) 2 ~ Vtt)) (2.2) 

where the scalar field <j) takes the role of a spectator field during inflation. Its potential is assumed to 
be of the form V{4>) = —3nMHj(f) and the coupling to the Maxwell tensor field, F^, is chosen to be 
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W{(j>) = e 2( ^^ M . Moreover, Hj = const, is the Hubble parameter and the scale factor is determined 
by a{rf) = —l/(j]Hi) for — oo < r] < nj < upto the end of inflation at 77/ . The resulting magnetic 
field is non helical and its power spectrum at the end of inflation on superhorizon scales is given 

by m 

Pb(M/)^ 1 * } (-pj HHim)\ (2-3) 

where 

n B =( 4 + 2 ^ +1 ) n < ^ (2.4) 

[ 4-2n n > 0. 

In the numerical solutions the magnetic field values are taken at present time. Thus assuming that 
the magnetic field only redshifts with expansion, so that Bi ~ 1/a 2 . Then the two point function 
of the magnetic field, 

(Bt(k,r,)B s {k!,rf)) = (<% " ^) (2-5) 

determines that the power spectrum scales as a -4 , so that at present time 

ciqPb(A;,?7o) = (Hirji)~ 4 PB(k,r}i). Then together with the parameters chosen as in [15], —r\i = 
1CP 24 Mpc for the end of inflation, ao = 1, and introducing a pivot scale ki, = 1 Mpc -1 for 
the magnetic field spectrum and a gaussian window function W(k, k m ) to model damping of the 
magnetic field on small scales thereby imposing an upper cut-off of the wave numbers, k m , the 
spectral function at present time is given by 

PM^ % ) = 10- 2 - 2 --ii^( I ^ T ) (A) W(k,k m )M V c-i (2.6) 



where W(k,k m ) = it 3 / 2 e \ k ™) , so that j d 3 kW(k,k m ) = kf^. The maximal wave number is 
given by [27] 

k m ~ 198.454 f Mpc" 1 , (2.7) 

for the values of the bestfit ACDM model of WMAP9, n b = 0.02264/i~ 2 and ^ = 0.7 [28]. Compar- 
ing the spectrum (|2.6p with the form of the magnetic field spectra used in, e.g., in [19] the spectral 
index there becomes ng — > rig + 3, so that the scale invariant cases correspond here to = and 
in [19] to tib* = —3. Thus the magnetic field strength today smoothed over the diffusion scale is 
determined by 

<p 2 (x,,o)> = io- 19 - 8 - 24 - 3 -r 2 r ) (j^-tJ 3 g 2 . (2.8) 



Defining the magnetic field strength B = y (B 2 (x, t]q)) and using equation (12 .7p yields to 

.^^(^ri^), ,, 9) 
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Figure 1: The root mean square magnetic field strength smoothed over the diffusion scale as a 
function of the spectral index ns- 



which is shown in figure [TJ Therefore for ub = 0.02 corresponding to ng, = —2.98 the smoothed 
magnetic field strength is of the order of 1 nG. The scale invariant case ub = corresponds to 
n = 2 or n = —3. As shown in [15J only the case n = 2 is allowed since in the other case strong 
back reactions would develop since the electromagnetic energy density dominates over the inflaton 
energy density. 



3 Cross correlation functions of primordial curvature and mag- 
netic fields 

In |15j the 3-point cross correlation between the magnetic field energy density proportional to S 2 and 
the fluctuations of the spectator field 8(j) has been calculated. However, since the anisotropics of the 
CMB are determined by the magnetic energy density contrast A# and the magnetic anisotropic 
stress ttb it is necessary to obtain the 3-point cross correlation functions allowing for arbitrary 
components of the magnetic field. Starting with the expression found for the cross correlation of 
the gauge potential A{ of the electromagnetic field and 5(f> [15] the 3-point cross correlation involving 
arbitrary components of the magnetic field is determined to be, in Fourier space, 



( m \ 1 ; Vl) B l (k2,m)B 3 (h,r 10 )) 
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where I\ and I 2 are integrals depending, in particular, on the parameter n which can be found in 
[15j . The special case of n = 2 will be presented in the notation used here below. This corresponds 
to a scale invariant magnetic field spectrum which is the cosmologically relevant case |15j . 

In order to relate the perturbations in the scalar field to the curvature perturbation the spectator 
field is identified with a curvaton field as suggested in |15j . In particular, we assume the simplest 
curvaton model [13] so that the total curvature perturbation £ is generated by the curvaton after 
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inflation and that any other contributions to the curvature perturbation are negligible. Moreover 
it is assumed that the curvaton dominates the energy density before its decay. Then following [13 

6</){k) 
M 



2 5)«*> 



(3.3) 



where the asterisk denotes the epoch when the mode k leaves the horizon during inflation assuming 
no nonlinear evolution on large scales [29]. Therefore calculating (C^Cp) a t the pivot scale of 



WMAP9 kp = 0.002 Mpc" 1 yields to 



(3.4) 



where the dimensionless power spectrum of two random variables X and Y, V^xy)^ defined by 

2n 2 



k 3 



(3.5) 



has been used. 

The magnetic energy density contrast and anisotropic stress are defined in terms of the photon 
energy density p 7 and pressure p 7 = |/3 7 . This is due to the weakness of the magnetic field which 
does not contribute to the background quantities (for a different approach see [30]). The energy 
density contrast A# in fe-space is defined by (for details, e.g., cf. [19] ) 



(3.6) 



So that one of the contributions to the cross correlation between the primordial curvature mode 
and the magnetic mode is determined by the two-point correlation function (£(/c')As(/c)). Taking 
the continuum limit Yl 

on superhorizon scales as given in [15] yields to 



,. -r J /g^ a and using the expressions for the integrals I\ and I2 for n = 2 
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The two-point correlation function {<^(k)irB(k')) is calculated using the expression for the 
anisotropic stress of the scalar mode (e.g., [19]) 



ir B (k) 



2p 7 o 



q i 
For the cosmologically interesting case n = 2 this yields to 
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The dimensionless power spectra determining the cross correlation between the curvature mode 
and the magnetic modes are shown in figure [2] (left) for S = 5 x 1CP 4 corresponding to the upper 
bound on this parameter inferred in [15]. As can be appreciated from figure [2] whereas the curvature 
perturbation is correlated with the magnetic energy density contrast, it is anticorrelated with the 
magnetic anisotropic stress. For comparison, in figure [2] (middle) the auto- and cross correlations 
of the magnetic modes are shown for magnetic field strength B = 1 nG and spectral index tib* = 
—2.98. The expressions for the correlation functions of the magnetic modes, V/a b A b ) (^)> ^V fl 7r s ) (^) 
and V(A B Tr B )(k) can be found, e.g., in [19] . Comparing the left and middle figures of figure [2] shows 
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Figure 2: Left: The dimensionless spectra determining the cross correlation function of the pri- 
mordial curvature perturbation and the magnetic energy density contrast as well as the magnetic 
anisotropic stress. Middle: The dimensionless spectra determining the auto- and cross correla- 
tions of the pure magnetic mode for B = 1 nG for a nearly scale invariant magnetic field, which 
corresponds to ub = 0.02 or ng» = —2.98 in the conventions used, e.g., in [19]. Right: The ratios 
determining the generalized Schwarz inequality for (C^s) and {(kb)- The numerical solutions are 
calculated for the parameters ^ = 5x 10 -4 , rji = — 10~ 24 Mpc [15] . 

that the amplitudes of the spectral functions of the cross correlations between curvature and the 
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magnetic modes are larger than the auto- and cross correlation functions of the magnetic modes. 
In figure [2] (right) the ratio 

Xs= r^C|L ; E = A b ,tt b (3.12) 

V^CCIIP<EE>I 

is reported which determines the generalized Schwarz inequality and as can be appreciated is well 
below the unity bound. 



4 Results 

The CMB angular power spectra are determined by the brightness function which in the line-of-sight 
approach [31] is written for each component, for the scalar mode, [32] 



®^ = jf* dns§ (k, r,)h [k( V0 - v)) (4.i) 

where Sq is the source function and X denotes £, A^j and 7r#. This determines the temperature 
auto correlation function 

C TT,(XY) =^y M ^jr % d^ <X y > (fc)^(fc,7,)^(fc^)(j / [fc( % -77)]) 2 (4-2) 

and similarly for the auto correlation function of the E-mode and the temperature polarization 
cross correlation. The resulting angular power spectra are calculated using a modified version 
of CMBEASY |33| . It is based on the modified version of [19J where the initial conditions and 
evolution equations including the magnetic field contributions can be found. In all solutions the 
best fit values of the 6-parameter ACDM model of WMAP9 [28] have been used, in particular, 
n b = 0.0463, U A = 0.721, V^kp) = 2.41 x 10~ 9 , n s = 0.972 and the reionization optical depth 
r = 0.089. Moreover, the nearly scale invariant magnetic field has field strength B = 1 nG 
and spectral index tib* = —2.98. The factor S determining the cross correlation functions is 
| = 5x 10~ 4 . 

In figure [3] the angular power spectra determining the temperature auto correlation of the CMB 
are shown. The contribution due to the cross correlation between the curvature mode and the 
magnetic modes dominates on large scales over the one due to the pure magnetic mode, that is 
the compensated mode for which the initial conditions are such that the contributions from the 
neutrino anisotropic stress and the magnetic anisotropic stress cancel each other. There is also 
the so called passive mode [20] due to the presence of a magnetic field before neutrino decoupling 
which causes an additional contribution to the curvature perturbation amplitude proportional to the 
magnetic anisotropic stress. However, since this is an adiabatic like mode it is not further studied. 
Similarly, the contribution due to the curvature-magnetic modes cross correlation dominates at low 
multipoles over the one due to the compensated magnetic mode for the polarization autocorrelation 
(cf. figured]) as well as the temperature polarization cross correlation (cf. figured]). In figure [6] the 
linear matter power spectrum is reported. On large scales the contribution due to the curvature- 
magnetic cross correlation dominates over that of the pure magnetic mode. On small scales the 
former contribution becomes subleading resulting in the characteristic rise of the linear matter 
power spectrum in the presence of a magnetic field. This is due to the presence of the Lorentz 
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Figure 3: Angular power spectra determining the temperature autocorrelation of the CMB due 
to the cross correlation between the curvature mode £ and the A# mode and due to the cross 
correlation between the curvature mode and the ttb mode (left). The contributions to the total 
angular power spectrum include the adiabatic mode (Adiab.) , the pure magnetic mode (PMM) 
and the total cross-correlated magnetic-curvature mode (CO) (right). These have been calculated 
for WMAP9 ACDM best fit parameters. 




Figure 4: Angular power spectra determining the autocorrelation of the polarization E-mode of 
the CMB due to the cross correlation between the curvature mode £ and the Ag mode and due 
to the cross correlation between the curvature mode and the ttb mode (left). The contributions to 
the total angular power spectrum include the adiabatic mode (Adiab.), the pure magnetic mode 
(PMM) and the total cross-correlated magnetic-curvature mode (CO) (right). These have been 
calculated for WMAP9 ACDM best fit parameters. 
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Figure 5: Angular power spectra determining the cross correlation of temperature and polarization 
E-mode of the CMB due to the cross correlation between the curvature mode £ and the A# 
mode and due to the cross correlation between the curvature mode and the ttb mode (left). The 
contributions to the total angular power spectrum include the adiabatic mode (Adiab.), the pure 
magnetic mode (PMM) and the total cross-correlated magnetic-curvature mode (CC) (right). These 
have been calculated for WMAP9 ACDM best fit parameters. 



term in the baryon velocity equation which on small scales dominates [2Q|, I19j . Thereby leading 
to a local maximum on small scales in the total linear matter power spectrum taking into account 
all contributions, that is the adiabatic and magnetic modes and their cross correlations. For 
wave numbers much larger than the maximal wave number k m determined by the damping of the 
magnetic field the linear matter power spectrum rapidly decays. In the case under consideration 
here k m = 198 Mpc" 1 (cf. equation flU]) ). 

5 Conclusions 

Observations of the CMB of unprecedented precision allow to test physics upto very early times. 
Observational evidence of large scale magnetic fields is abundant. There is however the open 
question of their origin. Assuming magnetic fields have been created in the very early universe, 
long before decoupling, they influence the spectra of the anisotropies and polarization of the CMB 
as well the linear matter power spectrum. There is clear observational evidence in the CMB pointing 
to the presence of an adiabatic mode which could be a primordial curvature mode generated during 
inflation (cf. e.g. [SI]). In general it has been assumed that the curvature mode and the magnetic 
field modes are uncorrelated. However, even at lowest order, that is not taking into account 
nonlinear effects due to the evolution of the magnetic field [22) . it seems natural to assume that 
non vanishing correlations exist depending on the particular model of magnetic field generation 
during inflation, such as the coupling of electrodynamics to the inflaton \24\ [25] or, as the case 
considered here, to a spectator field which is identified with the curvaton [15] . In [15] the bispectrum 
determining the cross correlation between the magnetic field and a scalar field was calculated. 
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Figure 6: The total linear matter power spectrum due to the pure magnetic modes (PMM), the cross 
correlation bwetween the curvature mode and the magnetic modes (CC) and the pure adiabatic 
mode (B = 0) (Adiab.) for the best-fit ACDM model of WMAP9. 

Identifying the scalar field with the curvaton here the cross correlations between the curvature and 
the magnetic modes is determined. Taking into account that the magnetic modes in the CMB are 
proportional to the magnetic field energy density and the magnetic anisotropic stress, respectively, 
both of which are quadratic in the magnetic field, the bispectrum or three-point function calculated 
in [15] induces a spectrum or two-point function relevant for the calculation of the CMB and matter 
power spectrum. The angular power spectra of the temperature anisotropies and polarization of 
the CMB have been calculated for a nearly scale invariant field with field strength of 1 nG. It has 
been found that the contribution due to the cross correlation between the curvature mode and the 
magnetic modes dominate at low multipoles over the contribution due to the pure magnetic mode. 
Similarly, the curvature magnetic mode cross correlated contribution dominates on large scales over 
the pure magnetic mode in the linear matter power spectrum. While on the contrary, on small 
scales the pure magnetic mode dominates, not only over the cross correlated contribution but also 
the adiabatic contribution which is due to the presence of the Lorentz term in the baryon velocity 
equation. 
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